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We present results from a numerical study of spherically- 
symmetric collapse of a self-gravitating, SU(2) gauge field. 
Two distinct critical solutions are observed at the threshold 
of black hole formation. In one case the critical solution is 
discretely self-similar and black holes of arbitrarily small mass 
can form. However, in the other instance the critical solution 
is the n = 1 static Bartnik-Mckinnon sphaleron, and black 
hole formation turns on at finite mass. The transition between 
these two scenarios is characterized by the superposition of 
both types of critical behaviour. 
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In a recent numerical study of gravitational collapse of 
a massless scalar field, a type of critical behaviour was 
found at the threshold of black hole formation [0 . More 
precisely, in the analysis of spherically-symmetric evo- 
lution of various one-parameter families of initial data 
describing imploding scalar waves, it was observed that 
there is generically a critical parameter value, p — p*, 
which signals the onset of black hole formation. In the 
subcritical regime, p < p*, all of the scalar field escapes 
to infinity leaving flat spacetime behind, while for super- 
critical evolutions, p > p* , black holes form with masses 
well-fit by a scaling law, Mbh oc (p—p*)''. Here, the crit- 
ical exponent, 7 ~ 0.37, is universal in the sense of being 
independent of the details of the initial data. Thus, the 
transition between no-black- hole/black- hole spacetimes 
may be viewed as a continuous phase transition with 
the black hole mass playing the role of order parame- 
ter. In the intermediate asymptotic regime (i.e. before a 
solution "decides" whether or not to form a black hole) 
near-critical evolutions approach a universal attractor, 
called the critical solution, which exhibits discrete self- 
similarity (echoing). Using the same basic technique of 
studying families which "interpolate" between no-black- 
hole and black-hole spacetimes, similar critical behaviour 
has been observed in several other models of gravitational 
collapse ||-|. 

In this letter we summarize results from numerical 
study of the evolution of a self-gravitating non-abelian 
gauge field modeled by the SU{2) Einstein- Yang-Mills 
(EYM) equations. In addition to its intrinsic physical 



interest, we have chosen this model since, in contrast to 
all previously studied models, it contains static solutions 
which we suspected could affect the qualitative picture 
of critical behaviour. Let us recall that these static solu- 
tions, discovered by Bartnik and Mckinnon (BK) [^,D, 
form a countable family X„ (n G N) of spherically- 
symmetric, asymptotically fiat, regular, but unstable, 
configurations. 
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FIG. 1. Schematic representation of "phase-space" for 
spherically symmetric Yang-Mills collapse, showing possi- 
ble end states of evolutions from a sufficiently general 
two-parameter family of initial conditions. The critical line 
00' demarks the threshold of black hole formation. An inter- 
polating family such as AA' exhibits Type I behaviour: the 
critical solution is the static BK solution Xi, and the small- 
est black hole formed has finite mass. Families such as BB' 
exhibit Type II behaviour: the critical solution is discretely 
self-similar (A « 0.74), black hole formation turns on at in- 
finitesimal mass, and mass-scaling with 7 « 0.20 is observed. 
At C, the two types of critical behaviour coexist. 

Our main new result is the observation that for certain 
families of initial data the static BK solution Xi plays the 
role of a critical solution separating collapse from disper- 
sal. Since in this case there is a finite gap in the spectrum 
of black hole masses, we call this "Type I" behaviour (in 
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analogy to a first order phase transition) , to distinguish it 
from "Type 11" behaviour (which we also observe) , where 
black hole formation turns on at infinitesimal mass. As 
is shown schematically in Figure 1, sufficiently general 
two-parameter families of initial data exhibit both types 
of critical behaviour. 

We consider spherically-symmetric Einstein- Yang- 
Mills equations with the gauge group SU{2). Follow- 
ing we write the general time-dependent, spherically- 
symmetric metric as 



(1) 



For the YM field we assume the purely magnetic ansatz 
which in the abelian gauge means that the field strength 
has the form |0] 

Wdt A O + W'dr A - (1 - W^)T3d^ A sin Mip, (2) 

where Q — Tid'd -\- T2sin'ddip, Ti (i = 1,2,3) are Pauli 
matrices, an overdot denotes d /dt and a prime denotes 
d/dr. Thus the matter content of the model is described 
by a single function, W{r,t)^ which we hereafter refer 
to as the Yang-Mills potential. We note that, as follows 
from (2), the vacua of the YM field are given hyW — ±1. 
We introduce auxiliary YM variables, $ = W and 11 = 
aW /a. The dynamics of the EYM model can then be 
computed using the following set of equations (see j?)): 
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W[T,i) = ±1 
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We have solved the initial value problem for many one- 
parameter families of asymptotically flat, regular initial 
data, some of which are listed in Table I. To ensure regu- 
larity at the origin we require that W{r,t) ±1 + 0{r'^) 
as r —^ 0. The numerical results described below were 
generated using a modified version of the adaptive-mesh 
algorithm used to perform the original scalar field cal- 
culations 0. The sensitivity of the mesh-reflnement al- 
gorithm was again very helpful in efficiently computing 
near-critical solutions, but was not nearly as crucial as it 
was for the scalar case. In fact, we have also reproduced 
most of the following results using a uniform-grid code 
much like the one described in 

As stated above, depending on the particular form of 
initial data used, we have found two types of critical be- 
haviour: 





FIG. 2. Late time profiles of marginally subcritical Type 
II collapse using family (a) (see Table 1). Note that the mass 
aspect, m{r,t), is defined via a — (1 - 2m/r)-\ The large 
number of echoes visible here (for fixed \p—p*\/p*), relative to 
the scalar case, is a reflection of the relatively small value of 
the echoing exponent (Aym ~ 0.74 versus Asf ~ 3.44). How- 
ever, this is partly offset by the fact that the mass-scaling ex- 
ponents for the two models also differ significantly. In general, 
dimensional/scaling considerations suggest that A5n = — 7^71 
where 5n and Stt are the changes in echo number, n, and 
TT = In |p — p*|, respectively [17]. 

Type II behaviour. — In this case we observe a contin- 
uous no-black- hole/black- hole phase transition, and the 
overall picture of criticality is very much analogous to 
that of scalar field collapse. For a generic Type II family, 
and in the near-critical, non-linear regime, we conjecture 
that the evolution asymptotes to a locally unique (up to 
r, t — !■ err, at, for arbitrary cr > 0) discretely self-similar 
solution, with an echoing exponent, A « 0.74. As with 
the scalar field case, we expect that the precisely critical 
solution echos an infinite number of times, exhibits un- 
bounded growth of curvature near r — 0, and is singular 
at the origin at some finite value of central proper time, 
T*. Figure ^ show profiles of various echoing quantities 
at r « T* from a family (a) calculation (see Table |) 
with \s — s*\/s* w 10~^^. Typical evidence for scale- 
periodicity is shown in Figure ^, using near-critical data 
from family (b) with parameter a varying. The level of 
agreement between the two independently computed es- 
timates of A provides us with a rough measure of the 
accuracy of our computation of the exponent. 

As expected, Type II families also exhibit mass-scaling 
in the super-critical regime: Mbh oc {p ~ p*)'^ . Typi- 
cal results are shown in Figure ^ and we estimate that 
the value 7 w 0.20 is accurate to a few percent. We 
note in passing that this result is another piece of the 
growing body of evidence which has shown that 7 is 
not constant across all collapse models [pp|jl0|]. As ar- 
gued in Ipdll, mass-scaling and universality (initial-data- 
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independence), strongly suggest that the stable mani- 
fold of a Type II solution, such as the one described 
here, is of codimension one. This picture, which pre- 
dicts that 7 is the reciprocal of the Lyapounov exponent 
of a single growing mode associated with a critical so- 
lution, has now been validated for at least two distinct 
matter models with continuously self-similar critical so- 
Perturbative treatment of criti- 



lutions lTT|jT|g,y,|i£ 

cal solutions with discrete symmetry is more involved, 
although considerable progress has been made for the 
scalar case Here we only remark that any tech- 

niques which work for the scalar field should be largely 
applicable to this solution. 
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FIG. 3. Illustration of scale-periodicity of Type II solution. 
This plot shows the superposition of a near-critical profile 
of $ (at a particular time) with the first two echoes which 
subsequently develop. The central proper time, To, at which 
the earliest profile is monitored is arbitrary; times Ti and T2 
and the rescaling exponent, Ap, are then chosen to minimize 
"I>(ln(r)-f nAp,T„) — <E>(ln(r),ro). An independent estimate of 
A is generated by first estimating the critical time, for the 
family, and then computing A.^ = ln((r* - T„)/(r* - T„+i)). 

Type I behaviour. — As stated previously, solutions in 
a Type I interpolating family asymptote to the static 
BK solution Xi as p p*. This behaviour appears 
generically for kink-type initial profiles of the YM po- 
tential W, such as family (c) in Table I. In this case it 
has already been established that the each of the BK so- 
lutions, Xn, has exactly n unstable modes (within the 
ansatz (2)); hence we are certain that the stable mani- 
fold of the critical solution has codimension one. Initial 
data with small \p~p*\ results in an evolution which ap- 
proaches Xi and stays in its vicinity for central proper 
time, T « —A In \p — p*\. The configuration then either 
disperses to infinity (p < p*) or collapses to a black 
hole with finite mass {p > p*). Typical results from 
a marginally subcritical Type I evolution are shown in 
Figure ||. 

Our calculations are somewhat complementary to 
those performed by Zhou and Straumann |^ . Those au- 
thors generally used initial conditions describing small 
deviations from the static solution Xi, and studied the 



subsequent evolution to verify the prediction of perturba- 
tive instability of Xi Here, by construction, we gen- 
erate the solution Xi as the boundary between collapse 
and dispersal, and thus immediately verify its instability. 
In this case, the reciprocal Lyapounov exponent of the 
single unstable mode yields the characteristic time scale, 
A, for the decay of Xi. Using central-proper-time nor- 
malization, the value computed from perturbation the- 
ory is A = 0.5519.... 
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FIG. 4. Mass-scaling of Type II solutions. Each marker 
type corresponds to a diflFerent family of super-critical com- 
putations. For each family, constants ai and /3i are chosen 
to unit-normalize the x-range and place the first data-point 
(smallest black hole) at the origin. For all plots, the least 
squares fit for the slope, 7, is 0.20, with an estimated un- 
certainty of a few percent. In addition, for all families, the 
unnormalized vr range was 18; thus in each case, the black 
hole mass spans a factor of e^^^ ~ 37. 

We can measure this exponent quite directly from our 
simulations by computing the variation of the lifetime 
of near-critical solutions with respect to variations in 
TT = \n\p — p*\. Specifically, defining Tr{n) to be the 
central proper time at which the zero-crossing of W(r, t) 
reaches radius r as it propagates outwards, we expect 
—dTr/d-K — > A as TT — ^ —00 and for sufficiently large r. 
Some numerical regularization is provided by monitoring 
Tr at several discrete radii ri^i = l...n and computing 
the averaged quantity TV(7r) = ^i-i (ti")- When this 

is done for family (c) with ri = 400, r„ = 475, n = 16, 
we find 0.5525 < -dfr/d-K < 0.5520 for -19 < tt < -10. 
In addition, we can get a good estimate of the unstable 
eigenmode by studying near-critical departures from the 
static solution. 

As noted in the introduction. Type I behaviour is 
clearly characterized by a gap in the black hole mass 
spectrum at threshold. We conjecture that the mass-gap 
is universal, and observe that our calculations suggest 
that it is very close (1% or less) to the total mass of the 
static solution (mi = 0.828640 . . .). 

Type I and II coexistence. — We can only briefly de- 
scribe what is one of the more interesting features of 
critical behaviour in the EYM model: the fact that for 
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certain two parameter families — such as (a) and (d) in 
Table | — there exists a critical line in parameter space 
(see Fig ||) which interpolates between Type I and Type 
II behaviour. The results we have obtained lead us to 
conjecture that the transition point (C in Fig |^) rep- 
resents coexistence of the two distinct critical solutions 
described above. In other words, near C, we see echo- 
ing occurring in context of the static, n = 1 background. 
Thus we can have have arbitrarily small black hole for- 
mation within a configuration which itself is arbitrarily 
close to forming a finite-mass black hole. We end by 
noting that the super-critical regime in this model, es- 
pecially in the cross-hatched overlap region sketched in 
Fig |l| is still very much terra incognita, although we have 
preliminary evidence that further phenomenological rich- 
ness lurks there. In particular, we have observed there 
an intriguing discontinuity in the spectrum of black-hole 
masses, which suggests the existence of another type of 
critical behaviour in the formation of black holes . 
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FIG. 5. Marginally subcritical Type I evolution. Here we 
plot the dynamical evolution of W{r,t) (solid line) and su- 
perimpose the static BM configuration Xi (dashed line). Ini- 
tially, the evolution (family (c), — « 10~^^) is nearly 
linear and almost purely ingoing. When the pulse arrives at 
the center, it sheds oH' YM radiation, approaches Xi and stays 
near it for some time, and then disperses to infinity. 
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TABLE I. Initial data families used in the calculations reported here. Listed are the family label, the form of the initial 
profile, W{r,0), family parameters, pi, nature of W's initial time derivative, and which behaviour is found by generating 
interpolating families. Initial data for family (a) is time-symmetric {W{r, 0) = 0); for all other families the initial configurations 
are almost purely ingoing. Finally, for family (a), the constants a and b are chosen so that VF(0, 0) — 1 and W^'(0, 0) — 0. 



Family W{r,0) pi W{r,0) Behaviour 

(a) (1 + o (1 + br/s) exp(-2(r/s)^)) tanh((x - r)/s) Ol 

(b) l + aexp(-((r-20)/5)«) a,S,q W'{r,0) II 

(c) (l-(r/<5)2)/((l-(r/<5)2)2+4r2)i/2 S W'ir,0) I 

(d) -l + 2aexp(-((r-17)/4)«) q,a W'lr,0) I, II 
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